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Thermodynamic analysis of an unconventional heat engine was performed. The engine studied has a 
number of advantages compared to state-of-the-art Stirling engines. The main advantage of the engine 
proposed is its simplicity. A power piston is integral with a displacer and a heat regenerator. It allows 
solving the problem of the high-temperature sealing of the piston and the displacer typical of all types 
of Stirling engines. In addition the design proposed provides ideal use of the displacer volume eliminating 
heat losses from outside gas circuit. Both strokes of the piston are working ones in contrary to any other 
types of piston engines. The engine can be considered as maintenance-free as it has no piston rings or any 
other rubbing components requiring lubrication. The only seal is contactless and wear free. It is located in 
the cold part of the cylinder. As a result the leakage rate in operation can be one-two orders of magnitude 
as small as that in Stirling engines. Balancing of the engine is much easy compared to Stirling engines 
with two reciprocating masses because of the only moving part inside the engine cylinder. The engine 
suits ideally to be fuelled with “difficult” fuels such as bio oil and can be used as a prime mover for 
micro-CHP systems. 

The thermodynamic model developed incorporates non-ideal features of the cycle, such as specific 
regenerator efficiency, dead volumes and other geometrical parameters of the engine. The model shows 
that the energy efficiency is highly sensitive to regenerator performance. For realistic geometric and 
operating parameters and the regenerator efficiency of about 95% the ultimate energy conversion effi¬ 
ciency of the engine proposed can be as high as 40-50%. 

A prototype of the engine was built and the feasibility of the engine concept was demonstrated. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

The strength of external combustion heat engines vis-a-vis 
internal combustion engines is compatibility with a wide variety 
of renewable energy and fuel sources. They may use a supply of 
heat from any sources such as biomass and biomass derived prod¬ 
ucts, municipal waste, nuclear, solar, and geothermal energy. Other 
important advantages of external combustion engines are low 
emissions due to continuous combustion and low noise, due to 
elimination of exhaust of high-pressure combustion products. 

Promising external combustion engine concepts are Stirling en¬ 
gines that convert thermal energy into mechanical energy of recip¬ 
rocating piston(s). The pistons are moved due to a cyclic change of 
gas phase working fluid pressure caused by its temperature and 
volume change. High thermal efficiency of the Stirling cycle (Car¬ 
not efficiency), long maintenance interval and fewer moving parts 
are additional advantages of Stirling engines [1], 
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In the operation of practical Stirling engines a significant devia¬ 
tion from the ideal cycle due to frictional losses, working fluid 
leakage, dead volumes, etc. takes place. Technical problems, in par¬ 
ticular balancing of pistons (or piston and displacer) reciprocating 
with a phase lag and a high-temperature sealing of the piston, cur¬ 
rently prevent the wide application of Stirling engines. 

Significant improvements of Stirling engines are expected by 
using engines in which the displacer is integrated with the power 
piston. In these engines the piston-displacer assembly is moved 
due to the cycling pressure change of the working fluid under 
the influence of change of its temperature and amount inside the 
working chamber. 

The purpose of this study was to perform thermodynamic anal¬ 
ysis of one of such engines. 

2. External combustion single-piston engine 

The engine studied in this article closely resembles the hot-air 
engine devised by Manson [2] and the machines patented by Bush 
[3] and Baumgardner et al. [4], Fig. 1 illustrates the basic principle 
the engine cycle. 
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Nomenclature 



c P 

working fluid heat capacity at constant pressure 

X 

displacement (position) of the displacer (m) 


(J/mol/K) 


distance between the displacer and the bottom at lower 

C v 

working fluid heat capacity at constant volume 


dead point (m) 


(J/mol/K) 

*2 

distance between the displacer and the top at upper 

/ 

frequency of piston reciprocation (s _1 ) 


dead point (m) 

L 

height of the displacer (m) 

x', 

=x,/h (-) 

M 

working fluid molecular mass (kg/kmol) 

*2 

=*2 lh (-) 

m 

working fluid mass (kg) 

a 

=T c /r„ (-) 

N 

number of moles [-] 

p 

= (S- 5 P )IS (-) 

N 0 

_ P 0 Sh 
— RT C ' ' 

y 

=£(-) 

H 

height of the cylinder (m) 

A = 

PI<x = PThITc A = PI a = pT H IT c 

h 

piston/displacer stroke (m) 



P 

pressure (Pa) 

Subscripts 

R 

gas constant (kJ/kmol/K) 

0 

ambient conditions (in the chamber) 

r 

regeneration efficiency (-) 

C 

cold, upper chamber 

T 

temperature (K) 

H 

hot, lower chamber 

S 

cross-sectional area of the cylinder/displacer (m 2 ) 

u 

upward motion 

s p 

cross-sectional area of the piston (m 2 ) 

d 

downward motion 


The engine consist of cylinder 1, divided by displacer 2 into hot 
(heated by combustion products) part la and cold (cooled by water 
jacket, for example) part lb. The displacer is not sealed; the annu¬ 
lar gap between the displacer and the cylinder provides a free gas 
flow between the hot and cold parts of the cylinder during recipro¬ 
cation of the displacer. The displacer is connected directly to power 
piston 3 so that they move together as one piece. The annular gap 
between the power piston and upper part of the cylinder is sealed 
by means of a contactless seal 4. The power piston is connected to 
the crank mechanism 5 to transform the reciprocating motion of 
the piston to shaft power. The cylinder is also equipped with two 
valves 6 and 7 shown in the upper and lower parts of the cylinder. 

The cylinder is filled with gaseous working fluid such as atmo¬ 
spheric air at ambient pressure. The engine starts when displacer 2 
is in the middle position (Fig. 1A). Heating of the lower part leads 
to increasing of air pressure inside the cylinder. Under the pressure 
acting on the power piston it moves upward driving crank mecha¬ 
nism 5. During this motion the displacer displaces air from the cold 
part of the cylinder to the hot one and the pressure gradually rises, so 
driving the piston further. When the piston reaches the end of the 
stroke (upper position, Fig. IB), valve 7 opens and hot pressurized 
air is released from the cylinder to atmosphere until the air pres¬ 
sures inside and outside the cylinder become equal. Then valve 7 
closes and piston is pushed down (Fig. 1C) by virtue of inertia of 
the crank gear. During this motion residual air is cooled and its pres¬ 
sure drops down well below the ambient pressure. As a result the 
outside atmospheric pressure accelerates the piston downward per¬ 
forming the vacuum or suction power stroke. At the end of the stroke 



Fig. 1 . Basic principle of the engine. 


(Fig. ID) upper valve 6 opens communicating the internal volume of 
the cylinder to atmosphere. Then valve 7 closes and the piston is 
pushed upward by the crank gear inertia displacing air in the cylin¬ 
der from the cold chamber to hot one. Pressure in the cylinder rises 
driving the piston upward and the cycle repeats itself. 

Instead of air any gas can be used as a working fluid. The engine 
can be modified to operate in a closed-cycle mode by using a con¬ 
tainer between the inlet and outlet valves. In such a closed cycle 
the pressure can be well above ambient pressure. This way the 
power generated can also significantly be increased. 

To increase efficiency a heat regenerator can be used. For not very 
high power ranges (1-5 kW) a well developed side surface of the dis¬ 
placer can be used as the regenerator. In that case the valves must be 
located only in the cold area of the cylinder (valves 6 and 7). 

3. Thermodynamic model and analysis of the heat process in 
the engine 

Fig. 2 shows a sketch of the engine together with the notations 
used for the thermodynamic analysis. 

3. J. Mass of the gas in the engine cylinder 

The mass of the gas inside the cylinder depends on the stage of 
the cycle, or on the direction of the piston/displacer motion. 

3.1.1. Upward motion 

The amount of the gas (mass or number of moles) in the cylin¬ 
der during the upward motion can be calculated at the lower dead 
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Fig. 2. View of the engine. 
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point when the gas pressure is known and equal to the ambient 
pressure or the pressure in the engine casing P 0 . 

Po = (,h + x 2 )Sp RTc 


W u = S(1 - p) [\p u (x) - Po) dx = S(1 - p) N u [" dx - P 0 h 
Jo Jo o - X 


= S(1 - P) 


N U A In 


B — h 


-Poh 




Nu = N uC + N uH = 


PoS 

R 


(x 2 + h)p x, 
T c + T h 


PoSh 

RT C 


((1 + x' 2 )p + x\a) 


= N 0 ((1 +x' 2 )p + x\ a) 
where N 0 =%*; x\=% ; x' 2 = x f. 

3A.2. Downward motion 

Accordingly, when the piston moves down the mass of the gas 
in the cylinder can be calculated at the upper dead point. 

Po = (x, + h)S RTH 


During the downward stroke the work done is 

W d = S(1 - P) [ (P d (x) - Po) dx = S(1 - P) L [ dx + P 0 h 
Jh L Jh O—X 


= S(1 - P) 


-N d A\n 


B 


+ Poh 


B-h 

The total useful work per cycle is 

W=W u + W d = S(l-/?)(N u -N d )Ahi^=S(l-/)ANAhi^ 
Using the expressions for AN, A, and B we get 

The work done can be rewritten as 

?-l 


fx'+l+x' 


W = P 0 hS( 1 - p) In 1 + 


1+^+x' 


N d = N d c + NdH = 


PoS \x 2 p 
R [t c 


Xi +h 

T h 


N 0 [x 2 P + {1 +xi)a] 


The difference in the gas amount in the cylinder during upward and 
downward movements is 


AN = N u - N d 


PpSh r p 
R [T c 


1 

Th_ 


= N 0 (P - a) 


One would expect the engine performs positive work if AN > 0, or 
ThP > Tc, or p > a. 

3.2. Gas pressure in the engine 


When the displacer is in the position x the gas pressure as a 
function of x can be expressed as 


P(x) 


Nh pT (N-N h ) 
S(x,+x) RTh Sp(x 2 + h-x) RTc 


( 1 ) 


where N is the amount of gas (number of moles) in the cylinder: 
N = N u - during the upward motion and N = N d - during the down¬ 
ward motion. 

From Eq. (1) we have 


Nh = — 

(x 2 + h 


NTr 


*>[& 


h 1 

(X 2 +ft-X)J 


Therefore 


P(x) = 


NRTrT, 


C‘H 


(^ + h-x)[^ + ^]s(x 1+ x) 

NRTcT h _ NA 

S[T H p{x 2 + h~x) + T c (x, + x)] ~ B - x 


It shows that the maximum mechanical work is achieved if 
Xj = x 2 = 0. In this case 

W,MS( (1-ffllnf 

If/ is the frequency of the reciprocation then the mechanical power 
of the engine is 

P w =jP 0 hS(l-/J)lnf (2) 


4. Energy efficiency 

The cycle is composed of four stages: 

1. The displacer moves up. The volume of the hot space 
increases from SX] to S(X] +/i). The amount of the gas in 
the hot space increases accordingly. 

2. The gas discharges from the cylinder. The piston/displacer 
is in the upper dead point. The volume of the hot space 
remains constant. The pressure drops to P 0 , the amount of 
the gas decreases. 

3. The displacer moves down. The volume of the hot space 
decreases from S(xt + h) to Sx^ The amount of the gas in 
the hot space decreases accordingly. 

4. The gas comes in/intakes to the cylinder. The piston/displac¬ 
er is in the lower dead point. The volume of the hot space 
remains constant. The pressure rises to P 0 , the amount of 
the gas increases. 

The general energy conservation equation and equation of the 
gas state are 

SQ T = dU + PdV ; T = ^ 


where 

RT c T h RT c h[P( \+x’ 2 ) + ax\] 

S[T H p - T c ] S[p- a]’ P-a 


where U = NC V T - internal energy; N - number of moles; Qj is the 
total amount of heat coming into the system (hot space). The tem¬ 
perature of the gas in the hot space is T H . 


3.3. Mechanical work performed by the piston and engine power 
During the upward stroke the work done is 


4.3. Stage 3 

The energy conservation equation for the gas in the hot space of 
the engine can be rewritten as 
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Qir - AN m C p T Un + Qi - C v T„AN m + S f P(x)dx 

Jo 


(3) 


0.2T = ~AN H 2 C p T H + 



dU = -A N H 2 C v T h 


where AN m the number of moles coming into the hot space, or the 
increase/change in the number of moles in the hot space during the 
1st stage; - the amount of heat transferred to the hot space dur¬ 
ing stage 11 Tlin - the temperature of the gas coming into the hot 
space. It depends on properties of the regenerator and heat ex¬ 
change between the gas and regenerator. 

In the ideal case the heat capacity of the regenerator (displacer) 
is much larger than the total heat capacity of the gas in the cylin¬ 
der, the heat exchange between the gas and the regenerator is 
instantaneous, the thermal conductivity of the regenerator along 
its axis is zero, and the temperature of the regenerator changes 
along its axis from T H to T c . 

Modeling of the heat exchange between the gas and displacer/ 
regenerator requires special consideration. Here we use simplified 
approach. 

The temperature of the gas coming into the hot space will be in 
between T c and T H 

Tun = IT Th + (1 — r\)Tc 

where r can be called the regeneration efficiency (0 ^ r < 1). 

From Eq. (3) and taking into account the calculations of the 
work above we get 

<2, = - A N m (C p T lin - C v T H )+SN u A\n B B _ h 

= -AN m ((C p r\ - C v )T h + C p (l - n)T c )) +SN U A in JL 


A N m = N m f -N mi 

where A N m is the change/increase in number of moles in the hot 
space during stage 1; N m f and N H u are the numbers of moles in 
the hot space in the end and in the beginning of stage 1. From the 
equations expressing the pressures in the hot and cold parts 


or 

0.2 = A N H 2 RT H 


a Nh 2 = N H 2i - N H2 f 

where A H H2 is the change/decrease in number of moles in the hot 
space during stage 2; N H2f and N H2i are the numbers of moles in 
the hot space in the end and in the beginning of stage 2. 

The decrease in the number of moles in the hot space during the 
gas discharge can be found from the equations: 

ThNh 2 i _ Tc(N u - N h 2 i) T H N H2 f _ Tc(Nd - Nmf) 

x, + h ~ x 2 p ' Xi + h ~ x 2 p 

As a result 


A N H2 = (N u - N d ) 


1 + x', 


= A N- 


1 +x' 


1 + x', + x' 2 A 1 + Xj + x' 2 A 


Therefore, the amount of the heat supplied to the hot space during 
stage 2 is 


AN(1 + x\)RT h 
1 + x', + x' 2 A ’ 


AN = N 0 (P - a) 


The regenerator will receive during stage 2 the heat 


0r2 = A N H 2 C P (T H - T 2out ) = A N H 2 C P (T H - (r 2 T c + (1 - r 2 )T„)) 
= AN H 2 C p r 2 (T H - T c ) 


where T 2ou , = r 2 T c + (1 - r 2 )T H is the temperature of the gas coming 
into the cold space.If the dead volume is neglected 

0 2 = Nq{P — ix)RTh 


0.2r = N o(/? — oc)C p r 2 (TH — Tc). 


ThNhu T c (N u - N mi ) T H N mf T C (N U -N mf ) 
x, (x 2 + h)/? ’ X1 + /1 X 2 p 


we find 


A N m = N u 


1 +xi 


1 + x\ + x 2 A (1 ~F x' 2 )A + x\ 


4.3. Stage 3 

The energy conservation equation for the gas in the hot space of 
the engine during downward movement of the displacer can be 
rewritten as 

r° 

0 3 i = ~A N H 3 C p T H + O 3 = -C v ThAN H3 +S P(x) dx 

Jh 


where 

A = p/a. 


03 = RT H AN H1 -SN d A\n B B _ h 


During stage 1 the regenerator will get the heat (basically take off 
the heat) 


On = -A N m C p (T Un - T c ) = -A N m C p u(T H - T c ) 
If the dead volume is neglected (x ^ = x 2 = 0) 

A N m =N U = = N 0 p 

K1 c 


and 

Qi = -N 0 p({C p ri - C v )T h + C p (l -r-i)T c )) +SN 0 pA In A 


Qri — ~NoPt^(T H — T c ). 


4.2. Stage 2 


A N H 3 = N h 3i - N H3 f 

where AN H3 is the decrease in number of moles in the hot space 
during stage 3; N H3i and N H3 f are the numbers of moles in the hot 
space in the beginning and in the end of stage 3. From the equations 

ThNh 3 j _ Tc(N,j - N H3 ,) T H N H3 f _ T c (N d - N H3 f) 

Xi " (x 2 + h)P 

_ A 

1 + Xi + Ax 2 yl(l + X 2 ) + Xi 

The regenerator will receive during stage 3 the heat 


Xi + h 

we find 

A M_M . 


X 2 p 


1 +xi 


Q r3 = A N H3 C P (T H - T 3om ) = A N H3 C P (T H - ( r 3 T c + (1 - r 3 )T H )) 
= AN H3 C p r 3 (T H - T c ) 


The energy conservation equation for the gas in the hot part 
during the gas release from the cylinder is 


where T 3ou , = r 3 T c + (1 - r 3 )T H is the temperature of the gas coming 
into the cold space. 
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If the dead volume is neglected 

A N H3 = N d = ^ = N 0 oc 

Kin 

and 

Q 3 = RT H N 0 a - ocSN 0 A In A 


Finally the thermodynamic efficiency of the engine is 

W 0-fling 

^(l-r)(l-a) + lnf 


(4) 


5. Discussion 


Qr3 = a NoC p r 3 (TH — Tc). 


4.4. Stage 4 


Q 4 t = AN H4 C p T 4in + Q 4 = J dU = AN/mCTh 

where T 4in = r 4 T H + (1 - r 4 )T c is the temperature of the gas coming 
into the hot space. 

0.4 = -AN H4 (C„T 4in - CJ H ) = -AN H4 ((C p r 4 - C V )T H + C p (l - r 4 )T c )) 


ANh4 = Nh4/ — Nh 4 i 

The increase in the number of moles in the hot space during the gas 
intake can be found from the equations: 


ThNh4i _ Tc(Nd - N H4i ) T H N H 4f 
x, " (x 2 + h)p ’ x, 

As a result 
AN H4 = (N U -N d ) 




x\ + (1 +x' 2 )A 


Tc(N u - N H4f ) 
(x 2 + h)P 


ANx\ 

x\ + (1 +x' 2 )A 


Therefore, the amount of the heat supplied to the hot space during 
stage 4 is 


‘ -*TfTT^M l(c ' r * - c ') r » + c »( 1 - 

The regenerator will receive (get off) during stage 4 the heat 


Q r4 = -ANn 4 C p (T 4 i„ - Tc) = —ANH 4 C p (r 4 Tn + (1 — r 4 )Tc — Tc) 
= —A N H4 C p r 4 (T H — T c ) 

If the dead volume is neglected 
Q 4 = 0 and Q r4 = 0 


The amount of gas in the hot space at the end of the cycle is the 
same as at the beginning of the cycle. Therefore 

AN m — A/V n2 — AN H3 + A/V,, 4 = 0 

The total amount of heat supplied to the regenerator is 

On = £ On = -A N m C p ri (T„ - T c ) + AN H2 C p r 2 (T H - T c ) 

i-1 

+ A N H3 C p r 3 (T H - T c ) - A N H4 C p r 4 (T H - T c ) = 

- C P (T„ - T c )(AN m ri - AN H2 t 2 ~ A N H3 r 3 + AN H4 r 4 ) 

Since = 0 we may assume that rj = r 2 = r 3 = r 4 = r. 

For simplicity we assume that the dead volume is zero. In this 
case the total amount of heat supplied to the hot space per cycle is 

i=4 

Or = E ^ = - N o/J((C p r - C V )T H + C p (l - r)T c )) + SN 0 (1A In A 

i=i 

+ N 0 {P - a)RT H + RThNqOC - aSNoAln/1 


or 

Qr = -N 0 /i((C p r-C„-R)r H + C p (l r)T c ))+SN 0 (P-x)A\nA = 
N 0 liC p (l-r)(T H -Tc)+SN 0 (P-<x)MnA 


The study performed shows that energy efficiency of the engine 
is determined by: (1) the properties of the working fluid; (2) the 
ratio of the diameters of the displacer and piston; (3) the ratio of 
the temperatures of the heater and cooler; (4) the efficiency of 
the heat regenerator. 

Figs. 3 and 4 show examples of the calculated thermodynamic 
efficiencies as a function of the ratio of the areas of the displacer 
S and piston S p . 

The model shows that thermal efficiency is highly sensitive to 
regenerator performance. In case of ideal regenerator (r= 1) the 
efficiency is 

r, 0 = 1 - j 8 = | 

Since the engine produces power if /? > oc the thermodynamic effi¬ 
ciency is always less than the ideal efficiency. It approaches the 



Fig. 3. Thermodynamic efficiency of the engine as a function of the geometrical 
parameter /? = (S - S p )/S for different regeneration efficiencies and a = T C /T H = 0.33. 



Fig. 4. Thermodynamic efficiency of the engine as a function of the geometrical 
parameter p = (S - S p )/S for different regeneration efficiencies and oc = T C IT H = 0.25. 
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Fig. 5. Mechanical power of the engine as a function of the geometrical parameter 
p = (S - S p )/S for different temperature ratios a. = T C IT H = 0.25 and 0.33. 

ideal efficiency if /i tends to a limit. However the engine power in 
this case tends to zero (Eq. (2)). 

Although the engine efficiency is always less than Carnot effi¬ 
ciency it can be as high as 50% for realistic geometric and operating 
parameters and the regenerator efficiency of about 95% Modeling 
of the heat regenerator and its optimization showed that its effi¬ 
ciency of about 95% can easily be achieved. Moreover, it can be 
made very close to 100%, say 98%. In that case the efficiency of 
the engine can be as high as 50-60% provided the geometrical 
parameter /? = (S — S p )/S has appropriate value. 

In view of its advantages energy losses in the engine are ex¬ 
pected to be much less than in Stirling engines. To estimate the 
influence of possible energy losses on the energy efficiency a de¬ 
tailed modeling of the engine was performed (to be reported). 
Due to difficulties with the calculations of the gas flow between 
the hot and cold parts of the engine the minimum temperature ra¬ 
tio a. = T c IT h was 0.54. The results showed that the difference be¬ 
tween the obtained thermodynamic efficiency, Eq. (4), and the 
expected practical efficiency is not dramatic. For example, at the 
temperature ratio of 0.54 the efficiency can be as high as 30%. 

Fig. 5 shows an example of the calculated power (Eq. (2)) for 
temperature ratios a = T C IT H = 0.25 and 0.33 as a function of the 


ratio of the areas of the displacer S and piston S p . The power was 
calculated for the following values of the engine parameters: inner 
diameter - 40 mm, piston-displacer stroke - 50 mm, pressure of 
the working fluid - 10 7 Pa, y = 1.667,/= 20 Hz. 

A prototype of the engine coupled with crank gear was built and 
the feasibility of the engine concept was demonstrated [5]. The in¬ 
ner diameter of the engine cylinder is 20 mm and the piston-dis¬ 
placer stroke is also 20 mm. The tests of the engine are under way. 

6. Conclusions 

The results show that optimal design and operating parameters 
should be selected as a compromise between power and efficiency. 
Very compact engine is capable of producing 2-4 kW power with 
practical efficiency above 30%. 

The engine studied has a number of advantages compared to 
state-of-the-art Stirling engines and suits ideally to be fuelled with 
“difficult” low-grade fuels such as bio-oil and can be used as a 
prime mover for micro-CHP systems. It seems to be an effective 
alternative and supplement to the established methods of energy 
conversion: internal combustion engines, microturbines, Stirling 
engines and fuel cells. 
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